Characterizing multiparticle entanglement in symmetric iV-qubit states 
via negativity of covariance matrices 
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We show that higher order inter-group covariances involving even number of qubits are neces- 
sarily positive semidefinite for N qubit separable states, which are completely symmetric under 
permutations of the qubits. This identification leads to a family of sufficient conditions of insepara- 
bility based on the negativity of 2k th order inter-group covariance matrices (2k < N) of symmetric 
iV-qubit systems. These conditions have a simple structure and detect entanglement in all even 
partitions of the symmetric multiqubit system. The observables involved are feasible experimental 
quantities and do not demand full state determination through quantum state tomography. 
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An important problem in quantum information the- 
ory [l|, 0] is the formulation of appropriate methods for 
detecting entanglement and then finding measures that 
quantify the degree of entanglement in multipartite sys- 
tems. These two issues are difficult to deal with in their 
full generality for examining multipartite systems and 
therefore, a strategy in their understanding is to focus 
on certain special symmetric states [3| . The choice of the 
states with specific symmetry is based both on feasible 
experimental possibilities and on mathematical consid- 
erations [J]. In this communication, we examine entan- 
glement properties of even number of qubits in quantum 
states obeying permutation symmetry. Symmetric mul- 
tiqubit states form an important class due to their ex- 
perimental significance [B|, Q . Taking advantage of the 
elegant mathematical structure associated with symmet- 
ric states, we propose a set of sufficient but not neces- 
sary conditions to detect entanglement via experimen- 
tally amenable inter-particle covariance matrix. The in- 
separability conditions obtained here are generalizations 
of our earlier result Q for pairwise entanglement in sym- 
metric multiqubit states. It is important to point out 
that in Ref. [7j these conditions are shown to exhibit a 
similar structure, involving the qubit cross-correlations, 
like those for the Gaussian states Q and thus our ap- 
proach reveals a structural parallelism between the con- 
tinuous variable states and multiqubits considered here. 

Symmetric multiqubit states remain invariant under 
any permutation of the qubits and are therefore re- 
stricted to a (N + 1) dimensional subspace of the en- 
tire 2 N dimensional Hilbert space - allowing for a sub- 
stantial reduction in the state space. This is the maxi- 
mal multiplicity space of collective angular momentum, 
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x,y,z, ( (Tpi is the Pauli op- 



Ji 

erator of /xth qubit) and is spanned by the eigen states 
{|f ,M) ; -f < M < f } of J 2 (with maximum 
eigenvalue J = N/2) and J z . 

An arbitrary N qubit system is characterized by the 



density matrix 

P = T^v J]] T ai a2 ... QJV (<7i ai (T2a 2 ■ ■ ■ ^Non) > (1) 

where u m = (7 (g> I ® . . . (g> a a (g) I ® . . .), - with a a 
appearing in the [i th position - denotes the Pauli oper- 
ator of the /x th qubit; ax, 0.2, ■ ■ ■ , ajv = 0, x,y, z and 
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(! . ! the real coefficients T ai a2 ... QJV are the 



averages 



= Tr [p (a lai a 2a2 ...erjvow)] 

= (ciaiO'2Q2 ■■■ a Na N ), 



(2) 



and X00...0 — 1 gives the normalization condition. Note 
that the total number of independent parameters in this 



A' 



-1 = 2^ 



1. For 



trace-class density matrix is 

multiqubit states obeying exchange symmetry, the state 
parameters T ai Q2 ... QJV are symmetric under interchange 
of any pair of indices (corresponding to swapping of the 
qubits). So, the total number of parameters reduce to 
(N + l) 2 — 1. Setting N — I indices equal to and re- 
maining / indices taking values x, y, z, we obtain moments 
of I th order (/ < N) : 



T, 
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1 ■ • • &Ui ) — ^-ii 12 ■■■H 00... (3) 

where ii, i 2 , ...,it=x, y, z. 



It is convenient to introduce collective multi-indices 

* = {h (2 ■■■ ik} , 3 = { Ji 32 ■■■ 3k }, 
so that the moments T( 2fe ) of even order 2k (with 
k = 1,2,..., [N/2]) may be arranged as 3 fe x 3 k real 
symmetric matrices and moments of k th order are 
arranged as 3 k componental columns: 
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Let us consider k qubit operators A*" ^ and B^, asso- 
ciated with two different groups a and b: 



A (k) - o 



B 



(fc) 



(T &iii cr &2. 



(5) 



Arranging A( fc ) and BW as a column £( fe ' of 2 • 3* opera- 
tors, (correspondingly, £}- k ^ — , B^ fc ^), as a row 
of operators), we define the 2fc th order variance matrix, 
as in Ref. Q 



V (2fc) = 1 ^^A^)+ h.c 



(6) 



where A£W = |( fc ) - (|^). Note that by construction 
©, y( 2fe ) is (2-3 fc x2-3 fe ) dimensional real symmetric pos- 
itive semi definite matrix. The elements of the variance 
matrix are 



V, 



(2k) _ 



1 



= ~ (uifuf })-{«d, «r>}j (7) 



( where £f } = ^ § J + Q ' Q l and V< 2fc ) is 

cast in a (3 fc x 3 fc ) block form, 

y(2fe) : 



^4(2fc) C(2fc) 
C (2fe)T g(2fe) 



(8) 



Clearly, the off-diagonal block C( 2fe ) corresponds to 2k th 
order covariances among the inter-group of multiqubits 
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(2k) 



rp{k) rp(k) 



(9) 



In the second line of we have used © and ([J|. 
The diagonal blocks A^ 2k ^ and are identical for 

a symmetric intra-group multiqubit system: A {2k) = 

(Af) because the intra- 
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group averages are the same viz., {A[ k) A (k) ) 
and (l? fe) ) = <S- fc) ) . Under identical local unitary trans- 
formations U ® U <g) . . . <S> U on the qubits - which preserve 
the symmetric space structure - the blocks of the variance 
matrix change as 



{B\ k) B {k) ) 



jS 2kS > -> A {2k) ' = 1ZA {2k) K T , 
c (2k)^ c (2ky = -jic {2k) ll T , 



(10) 



with iz = 



R® R< 



, a 3 fe x 3 k real orthogonal ma- 



k times 

trix, comprised of direct products (containing k factors) 
of 3 dimensional rotations R G SO (3) - corresponding 
uniquely to 2 x 2 unitary matrices U € SU(2). 

We now focus on the question: How would multiqubit 
entanglement manifest itself under different partitioning 
of a symmetric system? Our identification here is that 
the inter-group covariance matrix £( 2fc ) holds a key to 



symmetric multiqubit entanglement, coming from vari- 
ous even partitioning of the system. It is worth men- 
tioning at this juncture the important difference between 
the recent paper of Korbicz et. al. Q from our present 
work. These authors have proposed necessary and suffi- 
cient conditions for entanglement, reflected through two 
and three qubit partitions of a symmetric multiqubit sys- 
tem. Strikingly, the two qubit result is shown Q to be 
captured by the off-diagonal block of the variance matrix. 
An important open problem, concerning the inseparabil- 
ity features hidden in all the even qubit reduced systems 
of a symmetric iV-qubit state, is what we are addressing 
here, by generalizing our approach outlined in Ref. [7j. 

First of all, we note that positivity of the variance ma- 
trix y( 2k ) demands that the diagonal blocks A^ 2k ^ be pos- 
itive semidcfinite. However, there arc no constraints of 
positivity on the off-diagonal block C( 2fe ) as such, though 
separable symmetric states carry a distinguishing feature: 

Theorem: For every separable symmetric multiqubit 
state, inter-group covariance matrices C( 2fe ) of various or- 
der 2k < N are necessarily positive semidefinite. 

Proof: Consider a separable symmetric state of 2k 
qubits, which is decomposable as a convex sum of direct 



products of k qubit density matrices p. 



(fc). 
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l; 0<p w <i. (if) 



In this state the inter-qubit averages (A^ Bj) are also 
separable: 

W 

w 

= ^T^HTfM, (12) 

w 

where we have denoted Tr(p^ A^) = (A-i )w an d 
used the fact (AS } w = (B^ } w and the notation 
(A\ k) ) w = T^ k) (w). It is also clear that 

(A? \e P = <A (fc) )sep = T^iw). (13) 

w 

The real quadratic form Q (2fc) = x T c <2k) x = 
J2i j c s (2fc) Xt Xj, with an arbitrary real 3 fc componental col- 
umn X, when evaluated in the separable state (|11|) gives 

Qilp=J2p™ (^'h^) 2 -(e^ (t/V)*)) 

w \ w / 

(14) 

which is necessarily a positive semidcfinite quantity, im- 
plying in turn that C (2fc) > 0; k = 1, 2, . . . , [N/2], in a 
separable symmetric multiqubit state pip . This proves 
our theorem. □ 
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The above theorem leads to sufficient conditions for 
entanglement associated with even number of qubits in 
a symmetric state: If the inter-group covariance matrix 
C( 2fc ' is negative, then the symmetric multiqubit state 
exhibits 2k-qubit entanglement for k — 1,2,... with 
2k < N. This leads to a hierarchy of inseparability 
conditions, which test entanglement in even partition- 
ing. For k = 1, the condition < has been shown 
in Ref. to be a direct consequence of Peres-Horodecki 
partial transpose criterion [10( on two-qubit partitions 
of a symmetric multiqubit state. Thus, negativity of the 
covariance matrix serves as both necessary and suffi- 
cient for pairwise entanglement in the symmetric N qubit 
system. Any test which confirms the negativity of the 
real symmetric 3 fc x 3 k covariance matrix C( 2fc ) is suffi- 
cient to assert the inseparability of the symmetric multi- 
qubit state. In order to establish the negativity of C^ 2k ' , 
the Sylvester criterion ll| may be used: Negative value 
assumed by any of the principal minors of a hermitian 
matrix implies that the matrix is not positive semidefi- 
nite. Thus a series of sufficient conditions for entangle- 
ment of 2k qubits could be extracted from negative prin- 
cipal minors (of various orders) [12J of the corresponding 
covariance matrix C^ 2k ' . This brings out inseparability 
constraints involving a few correlation observables, mak- 
ing our criterion experimentally amenable. It may be 
noted that a series of inseparability conditions, resulting 
from negative principal minors of various orders, demon- 



strate 13j negativity of the (infinite dimensional) par- 
tial transpose of a bipartite continuous variable density 
matrix, which is the Peres-Horodecki criterion [l(| for 
infinite dimensional states. 

We now test our inseparability conditions C^ 2k ^ < 
by considering some well known examples of symmetric 
./V-qubit states such as, GHZ and W type states, which 
have attracted experimental focus [B| . 

For an even Q N qubit GHZ state [H|]: 

|GHZ N ) = -i= (|0jv> + |1jv>) = -j= (|00 . . .0} + |1 1 .. . 1» (15) 

we find that CW has one negative eigenvalue, 
A^~) = — 2(~~ 1 ) . The lower order covariances C^ 2k \ 
for k < N/2, are all positive semidefinite. This is obvi- 
ous because GHZ state is separable with the disposal of 
qubits. Thus, our C-matrix criterion is in concordance 
with the known result that the GHZ state is N -party 
entangled and is fragile under disposal of particles [151 ] . 

From experimental point of view, it may be noted that 
the lowest order (see [12]) principal minor, which records 
negativity of C'^ is the diagonal element, with the index 
i = {xxx . . . xy}: 



C)r> = t 



r-(^) 2 ={-: 3 



N/2 — even integer, 
N/2 = odd integer. 



(16) 



More specifically, even-A" qubit entanglement in GHZ 
states is revealed 16| by the measurement of the iV-qubit 



observable ( ) (where the 

qubit indices may be conveniently interchanged). 
Next, consider N qubit W-state (l5| : 



|Wjv) 



(|100 . . .0) + |0 1 . . .0) + . . .). 



(17) 



Here, the covariance matrices C^ 2k \ of all orders 
k = 1,2,..., [N/2], are negative ( with only one 
negative eigenvalue, A' - ) = — !^ y Therefore, W- 
state of N -qubits is confirmed to exhibit 2k qubit entan- 
glement for all values of k (with, of course, 2k < N). 
Here again, the 2k qubit entanglement is seen explicitly 
through the measurement of one of the diagonal elements 

. . . z}, for 



of the covariance matrix c| z 2fe ^ , with i = { 
which Cjf = T^ 2fc ^ = — 1. It is therefore sufficient to 
check that (<ri z o~i z ■ ■ ■ &2k z) is negative. Thus the W 
state has 2k qubit entanglement in all the even parti- 
tions 2k — 2,4, 6, . . . of the state. Our results confirm 
that the W state is robust under disposal of qubits [H| . 

We now investigate the implications of our insepara- 
bility conditions C( 2fc ) < for mixed states: To this 
end, suppose that experimentally produced W and GHZ 
states have noise-like admixture of incoherently super- 
posed symmetric states: 



(i 



where p n 



N+ 1 



P N + x |V) {ip\ , < x < 1, 



(18) 
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N 



m\ denotes the pro- 



jection operator onto the symmetric subspace 
Sym (C* 2 ® C 2 <g> ... ® C 2 ) of N qubits (P N 
is an identity matrix in the symmetric subspace of 
qubits and hence p n /(n + l) corresponds to a maximally 
disordered separable symmetric state), and \ip) is either 
a 7V-qubit |GHZ)at or |W)jv state. For the least eigen- 
value of the covariance matrix to be negative, the 
mixing parameter x has to be greater than a certain 
threshold value. We find the following range of x for 
which inseparability is indicated via neg ativity of C*W 
for N=2, 4 and 6 qubits: 

GHZ - noisy state : 0.25 < x < 1, for N = 2. 

0.0625 < x < 1, for jV = 4. 

0.014 < x < 1, for N = 6. (19) 

W - noisy state : 0.25 < x < 1, for N = 2. 

0.0899 < x < 1, for N = 4. 

0.042 < x < 1, for N = 6. (20) 

We observe that the x-range for iV-qubit entanglement 
is smaller for the noisy W state (see (|2TJ)) ). than that (fT^|) 
for the noisy GHZ state. But eventually for large N, 
both the noisy states remain entangled for all values of 
x. A more general trend (but a restricted domain for x) is 
found by examining the lowest order principal minor: the 
noisy GHZ state is N (even) qubit entangled, when < 
x < l, ( verified by demanding that the diagonal element 
; 0; the index i = {x x ... y}). For mixed 



noisy state of W, the inseparability range - for N qubit 
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entanglement - is identified to be (w + 2) < x < 1, resulting 
from the negative diagonal element t/™ ) = -x, with 
i = {z z ... z}. 

Entanglement in various even partitions of the W-noisy 
state (|T8ll are examined by using the n qubit reduced W 
noisy state, 



(N -n + l) 



Pn- 



+ x 



|Wjv_„) <W N -„| + — |0jv-,«) (0 N - n \ 



(21) 



The covariance matrices of all even partitions of the W 
noisy state are found to be negative, in a specific insep- 
arability range of the mixing parameter x. For example 
we find that a W noisy state is two-qubit entangled when 
N "* 12 < x < 1. Note that this inseparability range for two- 
qubit entanglement is much restricted than the one real- 
ized for entanglement in the largest even partition of the 
state (see (|20|) ). As N increases x — > 1 indicating that 
in the large N limit the two qubit partition of a noisy 
W state is separable throughout the range < x < 1. 
The n qubit reduced GHZ noisy state is a convex sum of 
three separable states p N - n /(N - n + l), |0jv- n ) (o w _„| and 
|ijv-n) (ljv-nl and is thus a separable state. 

In conclusion, we have here generalized our formalism 
of the symmetric two-qubit inseparability condition, ex- 
pressed in terms of inter-qubit covariance matrix [jj, to 
all even qubit partitions of symmetric ./V-qubit systems. 
This takes the form of a hierarchy of inseparability con- 
ditions on the inter-group covariance matrices of even 
order: < , k = 1, 2, . . . with 2k < N. Only 

for k — 1 (i.e., for two qubit partitions) the insepara- 
bility condition is both necessary and sufficient, and for 
all other values of k, these conditions are only sufficient. 
We have illustrated their use for both pure and mixed 
states involving GHZ and W type states. The symmet- 
ric multiqubit system considered here facilitates a richer 
analysis in terms of SO (3) irreducible tensors [T7[. The ir- 
reducible tensor approach leads to a family of criteria [l7| 
for entanglement based on covariance matrices involving 
collective angular momentum variables and is suitable to 
test inseparability in macroscopic atomic ensembles @. 
Our approach suggests further generalization to d-level 
symmetric multiparticle systems also. 
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whenever T^f^ 

assumed by any diagonal element I^ 2fc ' is 
fore sufficient to check negativity of 
For example, in an even iV qubit GHZ 
l£i = (o-i x o- 2x ...aN y o-N +lx o-N +2x ...o- Ny ) 
(where the multi-index % = {xx ... y}), while the 
A/2-qubit correlation observable {ai x c2x ■ ■ ■ o~ n_ y ) takes 
the value 0, 1 depending on whether A/2 = even or 
odd respectively. Measurement of N qubit correlation 
l u = (o-i x o-2 x . . .aN y ON +lx ON +2x . . .a Ny ) is 
therefore sufficient to verify multiqubit entanglement. 
A. R. Usha Devi, R. Prabhu, and A. K. Rajagopal, 
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